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The mixed spin-(l/2, 5*) Ising model on the union jack (centered square) lattice is investigated by establish- 
ing the mapping relationship with its corresponding eight- vertex model. An interplay between the nearest- 
neighbour interaction, the competing next-nearest-neighbour interaction and the single-ion anisotropy gives 
rise to a rather complex critical behaviour displayed in the reentrant phase transitions, the weak universal 
critical behaviour, as well as, a presence of first- and second-order phase transitions. The most interesting 
finding to emerge from the present study relates to a variation of the weak-universal critical exponents along 
the line of bicritical points, which is being twice as large for the mixed spin-(l/2, 5*) systems with the integer 
spin-S atoms as for the ones with the half-odd-integer spin-5 1 atoms. 
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1 Introduction 

One of the most fascinating exactly solvable lattice-statistical models exhibiting cooperative spontaneous 
order and extraordinary rich critical behaviour is being the eight-vertex model. Over the last three decades 
this model turned out to have various applications in seemingly diverse research areas, actually, its nu- 
merous special cases proved their usefulness as the ice models, the ferroelectric KDP and F models, the 
free-fermion models, the dimer models, or the familiar planar Ising models originally introduced for in- 
sulating magnetic materials |T). It is worthwhile to mention, moreover, that Baxter's exact solution of the 
symmetric (zero-field) eight- vertex model 1 2 1 has led to a remarkable breakthrough in the phase transition 
theory, since it provided a convincing evidence for a continuous change of critical exponents evidently 
contradicting the universality hypothesis 0. Exact results reported on the continuously varying critical 
exponents consecutively inspired Suzuki in order to propose the weak universality hypothesis |4|, which 
allows changes of the critical exponents that do not violate the concept based on scaling laws |5 1. 

Among the most famous special cases of the symmetric (zero-field) eight-vertex model exhibiting the 
non-universal critical behaviour belongs the spin- 1/2 Ising model on two square lattices coupled together 
by the four-spin interaction 0. As a matter of fact, this system shows the continuously varying critical 
exponents depending merely on the strength of the four-spin interaction [7|. Another interesting system, 
in which the non-universal critical behaviour has precisely been confirmed, is being the mixed spin-(l/2, 
S) Ising model on the union jack (centered square) lattice |8 9|. This mixed-spin system can be examined 
following the approach originally suggested by Lipowski and Horiguchi [8|, who solved one particular 
case of the mixed spin-(l/2, 1) system within the transformation to the corresponding eight-vertex model. 
Notice that the aforementioned procedure has recently been adapted by one of the present authors to explore 
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Oct = -1/2,+1/2 • S=-S,-S+l,...,S 

Fig. 1 Schematic representation of the mixed spin-1/2 and spin-5 Ising model on the union jack (centered square) 
lattice. Solid (broken) lines depict the nearest-neighbour (next-nearest-neighbour) interactions. 



another particular case of the mixed spin-(l/2, 3/2) union jack lattice J9). Despite the similarity of phase 
diagrams of both these systems, it is somewhat surprising that the mixed spin-(l/2, 1) model exhibits 
quite different variations of the critical exponents in comparison with its analogous mixed spin-(l/2, 3/2) 
version. With the aim of cumulating a set of exact results that could unambiguously clarify how the critical 
exponents may depend on the quantum spin number, the main purpose of present work is to provide further 
extension of the approach used previously in order to bring a deeper insight into the critical behaviour of a 
general mixed spin-(l/2, S) Ising model on the union jack lattice with arbitrary quantum spin number S. 

This paper is organized as follows. In Section|2] we provide a detailed formulation of the model system 
and subsequently, the main points of transformation method that ensures an equivalence with the eight- 
vertex model will be briefly mentioned. The most interesting numerical results for a critical behaviour will 
be presented and particularly discussed in Section[3]for two particular cases of the mixed spin-(l/2, 2) and 
spin-(l/2, 5/2) systems, respectively. Finally, some concluding remarks are drawn in Section^] 

2 Model system and its solution 

Consider the mixed spin-(l/2, S) Ising model on the union jack (centered square) lattice C schematically 
illustrated in Fig. [2 The mixed-spin union jack lattice consists of two interpenetrating sub-lattices A and 
B, which are formed by the spin-1/2 and spin-S sites denoted by empty and filled circles, respectively. The 
total Hamiltonian defined upon the aforedescribed lattice C reads: 

4N 2N N 

T~tmix = —J ^ S m (Ji — j' 2_j U kOl — S^, (1) 

(m,i)cj (fe,0cK m=l 

where cr, = ±1/2 and S m = —S, —S + 1, +S are Ising spin variables placed on the eight- and four- 
coordinated sites, J denotes the exchange interaction between nearest-neighbouring A — B spin pairs and 
J' labels the interaction between the A — A spin pairs that are next-nearest-neighbours on the union jack 
lattice C. Finally, the parameter D measures a strength of the uniaxial single-ion anisotropy acting on the 
spin-S sites and N denotes the total number of the spin-1/2 sites. 

In order to proceed further with calculation, the central spin-S 1 atoms should be firstly decimated from 
the faces of sub-lattice A. After the decimation, i.e. after performing a summation over spin states of all 
spin-S 1 sites, the partition function of the mixed-spin union jack lattice C can be rewritten as: 

u{<Ji,<Jj,<Jk,(Ji), (2) 

{a} i,j,k,l 
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where the summation is performed over all possible spin configurations available on the sub-lattice A and 
the product runs over all N faces (square plaquettes) each composed of one central spin-S* site surrounded 
by four spin-1/2 variables er^, aj, at, &i as arranged in Fig. [2 The Boltzmann factor u>(a, b, c, d) assigned 
to those faces can be defined as: 

+s 

cj(a,b,c,d)=exp[pj'(ab+bc + cd + da)/2] ^ exp(/?L»n 2 ) cosh[/3Jn(a + b + c + d)], (3) 

n=-S 

where (3 — l/(k-g,T), k-g, is Boltzmann's constant and T stands for the absolute temperature. At this 
stage, the model under investigation can be rather straightforwardly mapped onto the eight-vertex model 
on a dual square lattice £-p, since the Boltzmann factor oj(a, b, c, d) is being invariant under the reversal 
of all four spin variables. Actually, there are maximally eight different spin arrangements having different 
energies (Boltzmann weights) and these can readily be related to the Boltzmann weights of the eight-vertex 
model on the dual square lattice. If, and only if, the adjacent spins are aligned opposite to each other, 
then solid lines are drawn on the edges of the dual lattice Ct>, otherwise they are drawn as broken lines. 
Diagrammatic representation of eight possible spin arrangements and their corresponding line coverings 
is shown in Fig.|2] It can easily be understood that each of the eight possible line coverings around each 
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Fig. 2 Eight possible spin configurations around each central spin-S site and their corresponding line coverings at the 
vertices of dual square lattice. 

vertex of the dual lattice corresponds to two spin configurations, one is being obtained from the other by 
reversing all side spins. Since there is even number of solid (broken) lines incident to each vertex of the 
dual lattice Cv, the model under consideration becomes equivalent with the eight- vertex model. With 
regard to this equivalence, the partition function of the mixed-spin Ising model on the union jack lattice 
can be expressed in terms of the partition function of the eight-vertex model on the square lattice: 

Zmix(T, J, J', D) = 2Z & - V {U\,U>2, ... ,u> & ). (4) 

The factor 2 in this equation comes from the two-to-one mapping between spin and vertex configurations 
(two spin configurations always correspond to one vertex configuration). 

The Boltzmann weights, which correspond to eight possible line coverings of the eight-vertex model 
(Fig. 0, can directly be calculated from equation 

+s 

lu 1 = exp(/3J'/2) ^ exp(/3£>n 2 )cosh(2/3Jn), 

+S 

uj 2 = exp(-/3J'/2) ^ exp(/3£>n 2 ), 

n=-S 

+S 

0J3 = L04 = exp(/3Z?n 2 ), 

n— — S 

+S 

0J5 = ll>6 = ojj = dtp, = exp(pDn 2 ) cosh(/3 Jn). (5) 

n=~S 
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Unfortunately, there does not exist general exact solution for the eight-vertex model with arbitrary Boltz- 
mann weights. However, if the weights (0 satisfy the so-called free-fermion condition: 

UJ\UJ2 + ^3^4 = ^5^6 + ^7^8, (6) 

the eight- vertex model then becomes exactly soluble as the free-fermion model treated several years ago 
by Fan and Wu 1 10 1. It can readily be proved that the free-fermion condition holds in our case just as 
D — > ±00, or T — > 00. According to this, the infinitely strong single-ion anisotropy leads to the familiar 
phase transitions from the standard Ising universality class, since our model system effectively reduces to 
the simple spin-1/2 Ising model on the union jack 11 II or the square lattice 1 12| solved many years ago. 
Within the manifold given by the constraint (|6}, the free-fermion model becomes critical as long as: 

u>\ + u>2 + ^3 + ^4 = 2max{wi, u>2, 0J3, 014} . (7) 

It is noteworthy, however, that the critical condition ([7J yields rather reliable estimate of the criticality 
within the free-fermion approximation 1 10 1 even if a non-validity of the free-fermion condition (|6) is simply 
ignored. 

The second branch of exact solution occurs just as the Boltzmann weights satisfy the condition of 
the so-called symmetric (zero-field) eight-vertex model 1 2 1 : 

UJi = Ul 2 , ^3 = W4, UJ5 = £*>6, UJ7 = (8) 

As we already have = lu^, lu^ = ujq, and — uj%, hence, the symmetric case is obtained by imposing 
the condition uj\ = u>2 only, or equivalently: 

, R1 ,, ZlLseMPDni) cosh(2/3 Jn) 

exp(-/3J ) = — s — — . (9) 

E«=-s exp(/3Z?n 2 ) 

According to Baxter's exact solution [2 1 the symmetric model becomes critical on the manifold ^ if: 

wi + uj 3 + u 5 + lu 7 = 2max{wi, w 3 , oj 5 , w 7 }. (10) 

It is easy to check that lo\ always represents in our case the greatest Boltzmann weight, thus, the condition 
determining the criticality can also be written in this equivalent form: 



[ exp{(3 c Dn 2 ) + 2 ^ exp{(3 c Dn 2 ) cosh(/3 c Jn) 

S n=-S 

+s +s 
\ exp(/3 c L>n 2 )J [ ^ exp(/3 c Dn 2 ) cosh(2/3 c Jn) , (II) 

n=-S n=-S 



where C = l/(k-oT c ) and T c denotes the critical temperature. It should be stressed, nevertheless, that the 
critical exponents (with exception of 6 and 77) describing a phase transition of the symmetric eight-vertex 
model depend on the function fi = 2 arctan(w5o;7 /luiuj^) 1 ^ 2 , in fact: 

a = a' = 2--, /9=^-, v = v'=^-, 7=^, S = 15, r, = - (12) 



3 Results and discussion 

Before proceeding to a discussion of the most interesting results obtained for critical properties of two 
particular mixed spin-(l/2, 2) and spin-(l/2, 5/2) systems on the union jack lattice, let us briefly com- 
ment on their ground-state phase diagrams. Solid lines displayed in Fig. |3]represent ground-state phase 
boundaries separating five different long-range ordered phases, which emerge in the ground state when 
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Fig. 3 Ground-state phase diagram of the mixed spin-(l/2, S) Ising model on the union jack lattice for: (a) S = 2, 
(b) S = 5/2. Dotted rectangles schematically illustrate typical spin configurations within the basic unit cell of each 
phase, broken line connecting triple points shows a projection of the critical line lilt into J' — D plane. 



J > 0. Spin order drawn in dotted rectangles shows a typical spin configuration within the basic unit 
cell of each phase. It is quite obvious from Fig. [3]that both investigated spin systems have rather similar 
ground-state phase diagrams. Actually, when moving along the vertical axis in the up-down direction, 
then the single-ion anisotropy parameter D/J forces the central spins to lower their spin state, while when 
moving along the horizontal axis in the right-left direction, then the sufficiently strong antiferromagnetic 
next-nearest-neighbour interaction J' alters the structure of the ground state due to a competition with the 
nearest-neighbour interaction J. It is noticeable that the central spin-S* atoms are free to flip within the 
phases IV and V, where the interaction energy between spin-5 and spin- 1/2 sites effectively cancels out 
as a result of the antiferromagnetic alignment between the latter spins. Accordingly, the phases IV and V 
exhibit a remarkable combination of the spin order on the sub-lattice A (spin- 1/2 sites) and the disorder 
that appears on the sub-lattice B (spin-S* sites). There is nevertheless one principal difference between the 
two investigated spin systems: the central spins are non-degenerate within the phase V whenever the mixed 
spin-(l/2, S) system consists of the integer-value spin-S" atoms as they reside their lowest 'non-magnetic' 
spin state S = 0. Furthermore, it is also quite illuminating to see from Fig. [3]that there exist direct phase 
transitions between the phases I-V and II-V in the range of moderate values of D/J and J' / J, which are 
accompanied by a change of the spin order at the sub-lattice A as well as B. 

Next, our particular interest will be devoted to a broken line that depicts in Fig. [3] a projection of the 
exact critical line il It into the J' — D plane. As it can be clearly seen, this projection crosses T = 
plane along the first-order transition lines between the phases I-V and II-V, respectively. Referring to the 
studies reported on previously, this remarkable line of critical points can be identified as a line of bicritical 
points that have the non-universal continuously varying critical exponents [8 9 1. It will be shown later on 
that there indeed meet two second-order phase transition lines with one first-order transition line at each 
bicritical point and altogether, the line of bicritical points bounds a coexistence surface between different 
ordered phases. In this respect, we may conclude that the special critical line consisting of bicritical points 
determines a location of phase transitions between the phases I-V and II-V, respectively. 

Now, let us construct a global finite-temperature phase diagram. For this purpose, Fig.|4]displays the 
critical temperature as a function of the ratio J' /J for several values of the single-ion anisotropy strengths 
D/J. Critical boundaries depicted as solid lines represent exact critical lines obtained from the free- 
fermion solution (0, which can be applied on the variety (|6j fulfilled in two limiting cases D / J — > ±oo. 
On the other hand, dotted critical lines show estimated critical temperatures calculated from the free- 
fermion approximation simply ignoring a non-validity of the free-fermion condition (|6j for any finite value 
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of D/J. Second branch of exact solution, which is related to the aforementioned critical line of the 
symmetric (zero-field) eight-vertex model (lilt on the variety (|9), is displayed as a round broken line. 
Broken lines terminating at circled points show for completeness the zero-field condition (|9} calculated 
for several values D/J, since the zero-field condition should be identical to the first-order transitions 
terminating at bicritical points. As a matter of fact, there are strong indications supporting this concept: 
the zero-field condition (|9} should always show a coexistence of different phases and also in our case it 
always starts from and ends up on points giving the coexistence of the phases I-V and II- V, respectively. 




By combining the results illustrated in Figs. [3 a) and@Ja) we may come to the following conclusions: if 
the central spins are occupied by spin-2 atoms, then a right (left) wing of the displayed critical boundaries 
corresponds to the phase I (IV) if D/J > 0, while it corresponds to the phase III (V) if D/J < —2. 
Within the range — 2 < D/ J < the right wing always corresponds to the phase V, but the left wing 
corresponds either to the phase I or II depending on whether D/J is greater or smaller than the boundary 
value D^/J = —2/3. The situation exactly at this boundary value is also depicted in Fig. 0a), three 
critical lines should terminate at the special bicritical point (filled circle) in which three different phases I, 
II and V coexist together. The area bounded by the line of bicritical points and the zero-field condition for 
-Db/ J consequently separates two different regions of coexistence: the area on the lhs (rhs) from the zero- 
field condition determines a coexistence of the phases I (II) and V, respectively. The first-order transition 
line between the phases I (II) and V is shown for one particular case D/J = —1/3 (D/J = —3/2). 
Another interesting finding follows from a detailed comparison of Figs. |5Ja) and|4ja): in addition to 
the finite-temperature bicritical point (shown as empty circle) there appears another bicritical point which 
is located precisely at zero temperature when D/J = is selected. Notice that two finite-temperature 
bicritical points may even occur when D/J is slightly above zero, however, upon strengthening of D/J 
they move, along the line of bicritical points, closer to each other until they coalesce at D max / J = 0.0267. 

Similar finite-temperature phase diagram as discussed formerly for the mixed spin-(l/2, 2) Ising model 
on the union jack lattice exhibits also its analogous mixed spin-(l/2, 5/2) version (compare Figs. @Ia) 
and (b)). For the sake of brevity, we shall therefore not repeat its adequate description here and we just 
quote the most essential differences. Of course, the most interesting region of criticality is now restricted 
to another range of single-ion anisotropy where the phase coexistence emerges (—1 < D/J < 0). It is 
noteworthy that the midpoint of this interval (D^/J = —0.5) now determines a special point in which the 
coexistence of three different phases (I, II and V) occurs. Further, the doubling of bicritical points at fixed 
value of D/ J is even more pronounced for the mixed spin-(l/2, 5/2) system, actually, the coalescence of 
bicritical points occurs just as D max /J = 0.0385 is approached. Apart from these trivial findings, the 
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most significant difference can be observed in the exact critical line calculated in the limit D / J — > — oo. 
On the one hand, the reduced critical temperature (A;bT c / J) of the mixed spin-(l/2, 2) system increases 
linearly with J' /J, what means, that the critical temperature is independent of J and depends merely on 
J'. It should be realized that both phases III and V stable in the limit D/J — > — oo possess the central 
spin-2 atoms in their 'non-magnetic' spin state 5 = and owing to this fact, the model system with integer 
spin-5 atoms effectively behaves as a simple spin-1/2 square lattice. A slope of the T c = f(J') dependence 
indeed gives constant values k^T c / J' — ±[21n(l + \/2)] _1 for the phase III (upper sign) and V (lower 
sign) in accord with an exact critical point of the spin-1/2 ferromagnetic (antiferromagnetic) Ising square 
lattice 1 1 2 1 . On the other hand, the critical temperature of the mixed spin-( 1/2, 5/2) system taken in the limit 
D/ J — > — oo is strongly reminiscent to the one at D/J — > oo. Namely, both the dependences resemble the 
same critical line of the spin-1/2 Ising model on the union jack lattice II II : the dependence calculated for 
D / J — > oo is just rescaled and shifted by a factor of five due to five-times higher spin value of the central 
spins. Finally, we should remark a feasible appearance of the reentrant phase transitions observable in both 
the investigated spin systems. It is quite apparent that the reentrance appears on account of the competition 
between the nearest- and next-nearest-neighbour interactions, in fact, both the higher-temperature ordered 
phases IV and V exhibit the coexistence of a partial order and disorder, which has been conjectured as a 
necessary condition for appearance of the reentrant transitions 1 131 . 

Last, let us take a closer look at the variations of critical exponents along the line of bicritical points. At 
first, we shall investigate the line of bicritical points for the mixed spin-(l/2, 2) system. Fig. |5]displays a 
projection of this critical line into the J' — T c plane (the dependence scaled to the left axis) and simulta- 
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Fig. 6 Dependence scaled to a left axis shows how 
the critical temperature changes with the ratio J'/J, the 
curve scaled with respect to a right axis depicts variations 



of D/J along this line. Dotted lines are guides for eyes, of the critical exponent a along this line. _ 
neously, its projection into the J — D plane wnich is scaled with respect to the right axis. Next, Figs. |6J 

Qand|S]show how the critical indices a, j3, and 7, respectively, change along this critical line. It is quite 
evident from these figures that the exponents /3 and 7 monotonically increase upon strengthening of J'/J, 
while the critical exponent a shows a monotonous decrease as J'/J increases. Another interesting fact to 
observe here is that the critical exponents of special bicritical point (filled circle), at which the coexistence 
of three phases I, II and V occurs, are close (but not exactly equal) to their universal values. Finally, it is 
also worthwhile to mention that the bicritical points of the mixed spin-(l/2, 2) system exhibit a rather large 
variations of critical exponents ranging: a E (— 1,1), /3 s (1/16,3/16) and 7 S (7/8,21/8). 

For comparison, Figs. 191 121 illustrate relevant dependences of the critical temperature and critical ex- 
ponents for the mixed spin-(l/2, 5/2) system. As it can be clearly seen, the critical exponents (3 and 7 
approach their smallest possible values 1/16 and 7/8 by reaching both triple points with zero critical tem- 
perature, while the critical index a approaches there its greatest possible value 1, With high certainty, the 
same critical exponents of both triple points can be related to a similar spin alignment of three coexisting 
phases I-II-V and III-IV-V. As a matter of fact, this conjecture is in a good accord with what is observed 
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in the mixed spin-(l/2, S) system with integer spin-S* atoms, where the critical indices at both triple points 
differ significantly due to a presence of 'non-magnetic' state S = energetically favored in the phases 
III and V. It is worthy to notice, moreover, that the greatest values of critical indices (3 and 7, as well as, 
the smallest value of critical exponent a, are being almost completely equal to the values predicted by the 
universality hypothesis for planar Ising systems. With all this in mind, we may conclude that the mixed 
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exponents as they are restricted to smaller ranges: a £ (0,1), /?E (1/16,1/8) and 7 S (7/8,7/4). 



4 Concluding Remarks 

The present article deals with a critical behaviour of the mixed spin-(l/2, S) Ising model on the union jack 
lattice that represents one of few exactly soluble planar Ising systems accounting also for the interaction be- 
yond nearest neighbours. It is worthwhile to remark that the relatively precise and concise information on 
criticality of the considered model system has been obtained through an establishment of the mapping rela- 
tionship with its corresponding eight-vertex model, which has exactly been solved in some particular cases 
when its Boltzmann weights satisfy either the free-fermion condition (|6}, or the zero-field condition l|8}. 
In agreement with our expectations, an interplay between the nearest-neighbour interaction, the compet- 
ing next-nearest-neighbour interaction and the single-ion anisotropy gives rise to a rather complex critical 
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behaviour displayed in the reentrant phase transitions, the weak universal critical behaviour, a presence of 
first- as well as second-order transitions and so on. 

The most challenging question concerning with the critical behaviour of the considered model system 
should be connected with the remarkable line of bicritical points, which bounds a coexistence surface 
between phases with different spin alignment. It is noteworthy that the results presented in the pioneering 
work of Lipowski and Horiguchi |8| imply two times greater variations of critical exponents along the 
line of bicritical points for the mixed spin-(l/2, 1) union jack lattice when comparing them with the results 
recently reported for its analogous mixed spin-( 1/2, 3/2) version 1 9 1 . Owing to this fact, the main aim of the 
present work was to clarify how the critical exponents may depend on the quantum spin number S. In view 
of the exact results presented here it is quite tempting to conjecture that the mixed spin-(l/2, S) models 
may exhibit different variations of critical exponents depending on whether spin-S 1 sites are occupied by 
the half-odd-integer spins, or the integer ones. Of these systems, the mixed-spin systems with the integer 
spin-S 1 atoms have been found to exhibit greater changes of the critical exponents. 
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